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1 Introduction 

There are many different mathematical approaches to classification and to investigation of the super- 
integrable systems with very extensive literature on the subject, see some recent rewiews [51 1121 [T5] . 
Here we want to discuss one of the oldest but almost completely forgotten approaches, which relates 
superintegrable systems with the addition theorems. 

The story began in 1761 when Euler proposed a generic construction of the additional algebraic 
integral of equations of motion 

Ki dXi K2 dx2 

where P is an arbitrary quartic and k's are integer [5]. In fact, the Laplace integral for the Kepler 
model and entries of the angular momentum for the oscillator are particular cases of the Euler integral. 

Remind, that with geometric point of view Euler (and later Abel in his famous theorem) studied 
the points of intersection of the hyperelliptic curve /i^ — P{x) with any arbitrary algebraical curve 
whose equation in a rational form may be written as F(xi, . . . ,Xn) = (see the following symbolic 
picture for the Euler two-dimensional case) [3] . 




In the Jacobi separation of variables method these points are the points of intersection of the sepa- 
rated relations pj = P(xj) defined on the whole phase space with the classical trajectory of motion 
F{xi, . . . , a;„) = on its configurational subspace. 

So, Euler proposed algebraic construction of the classical trajectories of motion without any inte- 
gration and inversion of the Abel map. This particular algebraic construction of the trajectories of 
motion was one of the starting points of the Jacobi investigations of the last multiplier theory and of 
the inversion of the Abel map. 
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Such as algebraic curve F{xi, . . . ,Xn) = has many different parameterizations Xk{t), there are 
many different superintegrable systems associated with this curve. For instance, the overwhelming ma- 
jority of the known two-dimensional superintegrable systems are related with the different parametriza- 
tions of the conic sections. 

Remark 1 For the ellipse 

F(xi,X2) = ^ + ^-l = 0, 
we can use well known parametrizations by trigonometric or elliptic functions 

Xi = asinf, X2 — boost or xi = asnt, X2 = bent 
which may be associated with the oscillator and the Kepler problem, respectively. 

1.1 Additional integrals of motion and angle variables 

A classical superintegrable system is an integrable n-dimensional Hamiltonian system with n + k single 
valued and functionally independent integrals of motion on the whole phase space M . At /c = n — 1 
we have maximally superintegrable systems with the closed trajectories of motion. 

The Liouville classical theorem on completely integrable Hamiltonian systems implies that almost 
all points of the manifold M are covered by a system of open toroidal domains with the action-angle 
coordinates / = (/i, . . . ,/fc) and uj = (wi, . . . ,w„): 

{Ij,Ik} = Wi,(^k} = 0, {Ii.i^j} = (1-1) 

The n integrals of motion Hi, . . . ,Hn are functions of the action variables Ji, . . . ,/„ such that the 
corresponding Jacobian does not equal to zero 

detJ^O, where J,, = ^^^{h I-a) ^^3) 
The n ~ \ additional integrals of motion 

-^J ^fe' {i?i>j} = 0, j = 2,...,n, (1.3) 

fc 

are functionally independent on Hi{I), . . . , Hn{I) and any completely integrable system is maximally 
superintegrable system in a neighborhood of every regular point of M [24] . 

In generic case the action variables Ik are single-valued functions, whereas the angle variables Uk 
are multi- valued functions on M. However, additional functionally independent integrals of motion Ki 
have to be some functions on (jjj , which are sums of the multivalued functions ujk ■ 

So, for any superintegrable system we have the following question: 

How are single-valued functions obtained from the sums of multi-valued func- 
tions? 

In fact for all the known superintegrable systems we have only one answer: 

We have to apply addition theorem and to get one multi-valued function on 
the single-valued argument, which is the desired integral of motion: 

k 

So, the addition theorems could help us to classify algebraically superintegrable systems and vice versa 
we can get addition theorems by using known superintegrable systems. 
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2 The action-angle variables for the open lattices. 



In 1975 Moser proposed construction of the action-angle variables for the open Toda lattice. This 
construction may be extended on other integrable systems related with the quadratic r-matrix algebras 

For all these systems the action variables Ij are zeroes of the polynomial 

n 

A{X)^r + Hi\"-^ + ---Hn = l[{X-Ij), {^(A),A(/i)} = 0, (2.1) 

j=i 

whose coefficients are independent integrals of motion Hi , . . . , Hn in the involution. 

The angle variables ujj are defined by the second polynomial B{\) having the standard Poisson 
bracket with the first polynomial A{fj,) [53] 

{BiX),Aip)} - ^ {B{X)A{^i) - Bi^,)AiX)) , 77 e C, 

or 

{BiX), A{^l)} = (B{X)A{^l) - B{^l)A{X)) + (b{X)A{p) + Bi-^,)A{X)) . 

In this case angle variables cuj are equal to 

c^, =7?-MnS(/,), j = l,...,n, 77 eC, (2.2) 

Usually brackets between these two polynomials are either part of the Sklyanin algebra or the reflection 
equation algebra. It allows us to use the representation theory in order to get single-valued function 
A{X) and B{X) on M in initial physical variables, see [24] . 

Using the angle variables cuj (|2.2p we can obtain additional algebraic integrals of motion 

K, ^ e^^' - B(/,), {4, if,} - 0, J, 

for these maximally superintegrable systems. 

Among such maximally superintegrable systems we can find the open Toda lattices associated with 
the classical root systems An, BCn, and I?„ , the XXX Heisenberg magnets with boundaries, and the 
discrete self-trapping model with boundaries [24] . 

Remark 2 The properties of the trajectories of motion is independent on the choice of integrals of 
motion in the involution /i, . . . , /„ or Hi, . . . , H„. This choice affects on the parametrization only, i.e. 
on the explicit dependence of the coordinates qj on the parameter t. 

Remark 3 The global action-angle variables, the corresponding Birkhoff coordinates and the KAM- 
theorem for the periodic Toda lattices are discussed in [9]. Up to now the computer experiments 
indicate that periodic Toda lattice can not be superintegrable system on the whole phase space. 



3 The action-angle variables for the Stackel systems. 

The system associated with the name of Stackel [THl [10] is a holonomic system on the phase space 
M = R^", with the canonical variables q — {qi, . . . , qn) and p = {pi, . . . ,Pn)'- 

n 

= ^ dpj A dqj , {Pj,qk} = Sjk ■ (3.1) 

The nondegenerate nx n Stackel matrix S, whose j column depends on the coordinate qj only, defines 
71 functionally independent integrals of motion 

n 

Hk=J2{S-'hk{p',+U,{q,)). (3.2) 
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From this definition one immediately gets the separated relations 

71 

p) = Y,HuSu,-UMj) (3.3) 



k=l 

and the angle variables 



i=i 



It allows reducing solution of the equations of motion to a problem in algebraic geometry [20] . Namely, 
let us suppose that there are functions fij and Xj on the canonical separated variables 

= Ujiqj)pj, \j = vj{qj) , {q^:PJ} = S,j, (3.4) 

which allows us to rewrite separated equations (j3.3p as equations defining the hyperelliptic curves 

C, : = P, (A, ) EE u2 (A^. ) 1^^^ HkSk, (A, ) - U, (A, )^ , (3.5) 

where Pj {Xj ) are polynomials on Xj . In this case the action variables Ik — p.2p have the canonical 
Poisson brackets p.ip with the angle variables 



j^l ^d^dX,=j2^.,{p,,q,), (3-6) 



which are usually the sums of integrals "dij of the first kind Abelian differentials on the hyperelliptic 



curves Cj (|3.5p 



Remark 4 If the separated relations are obtained from the hyperelliptic curves by some substitution 
(13. 4p the Lagrangian submanifold of the corresponding Stackel system is a direct product of these 
hyperelliptic curves 

= Ci X C2 X • • • X C„. 

If these curves are equal Cj = C then JF may be identified with the Jacobian J(C) of C [20] . 
3.1 Addition theorems 

In generic case the action variables (|3.6p are the sum of the multi- valued functions iSij . However, if we 
are able to apply some addition theorem to the calculation of uji (|3.6p 



Ui = ^ ■dijipj,qj) &i{Ki) + const, (3.7) 

where Qi is a multi- valued function on the algebraic argument Ki{p, q), then one could get additional 
single- valued integral of motion Ki . 

Treating A„+i as a function of the independent variables Ai, . . . , A„ the addition theorem (j3.7p may 
be rewritten in the integral form 



y/P~{X) J VPn+lW 



or in the differential form 

'S'ii(Ai)dAi S'i„(A„)dA„ S'(A„+i) dA„+i /o o^ 

£1 = H • • • + Cn = £?i+l — U, £i — 

\/Pl{Xi) \' Pn{Xn) V-fri+l(A„+l) 

where choice of the is related with the choice of the contours of integration. In this case additional 
integrals of motion are X„+i and s{Xn+i) ^ which have to be single- valued functions on the whole 

V -Pn + l(An + l) 

phase space. 
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Remark 5 If we add the necessary words to the equation (|3.8p at Pi = Pj, for all the i,j — 1, . 
we obtain the famous Abel theorem, see as an example P^ . Recall, that Abel reduced equation 
to an algebraic identity (trajectory of motion), which is a consequence of an expansion into partial 
fractions. 

Another standard form of the addition theorem for the function /(x) looks like 

F[f{x), f{v)J{x + y)) - , Vx, y e M, (3.9) 

where F is some function associated with /. For instance, we can show different forms (|3.7p and p.9p 
for the simplest addition theorem 

ln(x) + ln(y) = In {xy) or e^e^ = e^+^ (3.10) 

associated with the zero-genus hyperelliptic curves Cj 

Cj : fi'j = ^'j(^j) = /j^j + .9j^j j = 1, . . . , n, (3.11) 

where fj, gj, hj are linear functions on n integrals of motion Hi, . . . , H„. 

Namely, if the separated relations are obtained from the zero-genus hyperelliptic curves (j3.1ip and 
Sij — Kij are constants we have 

n n „ 

i=i i=i-^ \Jfj^j+9j^3+hj 



i=i 



„ , 2/, A, +5, 



So, the angle variables uji are logarithms as for the Moser systems (|2.2p . Of course, if we want to get 
algebraic or polynomial integrals of motion Ki we have to impose some constrains on and fj. 

We have to underline that we use different hyperelliptic curves Cj in this example in contrast with 
the Euler example. 

3.2 Trajectories of motion 

According to the standard Jacobi procedure of inversion of the Abel map the variables qj are found 
from the equations 



where Pi — t and = aj. From these equations one gets trajectories in parametric form qj{t), i.e. 
coordinates qj are functions on the time Pi — t and parameters ai, . . . , an, P2 ■ ■ ■ , Pn- 

On the other hand, for the superintegrable systems if we substitute momenta from the separated 
equations ()3.3p 



Pj = 



J^'^kSkjiqj) - Uj{qj) 
\ k=l 



into additional algebraic integrals of motion Kj one gets trajectories of the motion in the algebraic 
form 

Fj ^ Kj{qi,. . . ,qn,ai,. .. ,an) - kj ^0, j = l,...,n-l. (3.12) 
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Remark 6 There are many different parameterizations of the algebraic curves. Since, trajectories 
may be associated with the different Stackel systems p.l2p . which are related by the canonical trans- 
formation of the time t ^ t: 

dt = v(q)dt, v(q)^'^^^. (3.13) 
detS" 



Here S and S are the Stackel matrices for the dual systems with common trajectories [2TJ[23] . Existence 
of the such dual systems is related with the fact that the Abel map is surjective and generically injective 
if n = g only. 

Of course, such dual Stackel systems have so-called equivalent metrics. 

4 Logarithmic angle variables 

In order to classify two-dimensional superintegrable systems associated with the addition theorem 
(j3.10p we can start with a pair of the Riemann surfaces 

C,: ^2^P,(A) = /A2+5,A + /ij, j = 1,2, (4.1) 

where 

f = aHi+ f]H2 + 7, gj = a'jHi + /3|i/2 + 7|, hj = a']Hi + P';H2 + 7^, 

and a, [3 and 7 are real or complex numbers. 

We can use addition theorem p.lOp if we we fix last row of the Stackel matrix only. Namely, 
substituting 

52,(A) = «:, (4.2) 

into the (|3.3ll3.6p one gets 

^ J S^A^Ml^ -J^- ^.r^'^ In (m. + ^^^) , (4.3) 



so the angle variable 

uj-2 = — = In 

Vf 



dp,{\) 



\=Vj {qj ) 



is the multi-valued function on the desired algebraic argument 



^1 \ ( „ „, I 2 



If Ki^2 are positive integer, then 

/ 1 \ ^i+^s 

(i^^ + V7^m) (4.4) 

is the generating function of polynomial integrals of motion Km and Ki of the m-th and m ± 1-th 
order in the momenta, respectively. As an example we have 

K,n = 2 {piUi P2 +P2U2 Pi), Ki = 1, K2 = 1, 

Km = 2P^ (2P2U2 Pi + PlUi P^^ -I- SpiUiplulf , Ki = l, K2 = 2, (4 5) 

Km = 2Pf (3p2U2Pi + PlUlP^j + Sp^ul (p2U2Pi -f ipiU^P^ /, = 1, ^3 = 3, 

where m = 1,3, m = 3,5 m — 3,7, because Pi 2 and / are linear functions on Hi^2, which are the 
second order polynomials on momenta. 
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The corresponding expressions for the Ki look like 
Ke = P[P!,^ + 4fp2U2{p2U2Pi + 2piUiP^) , 



Ke = P[P^' + l2P!,p2U2{P!,PiUi + P[p2U2)f + l%piUiplulP 



Kl = 1, K2 = 1, 
Kl 1, K2 = 2, 
Kl = 1, K2 = 3. 



(4.6) 



Of course, we can try to get another polynomial integrals of motion using the following recurrence 
relations 



K,_^^{K,,H2},, 



and 



{iJi, ifj} = , j = m, m ± 1, m ± 2, . . . . 



(4.7) 



see examples in pS . 

Remind, that we get polynomial integrals of motion Km and K( for the special Stackel matrices 
only. The imposed necessary condition (|4.2p leads to some restrictions on the substitutions 

In fact, after these substitutions we obtain the following expression for the Stackel matrix 
/ avf + afvi + a'l au| + W2 + 0^2 \ 

S{q) = 



"1 "2 



V 



det 5* 7^ 0. 



(4.8) 



So, for a given ki^2 expressions for t?2j (14. 3p yield two differential equations on functions u,v and 
parameters /3: 



S2j{qj 



Uj{q.j) 



KjUjv'j = + + Pf, j = 1, 2. 



(4.9) 



There are equations on u, v and parameters f3 because we have to solve these equations in the fixed 
functional space [25] , 

Proposition 1 // kj ^ equations |^.g| ) have the following three monomial solutions 

3! 



I 

II 
III 



/3 = 0, /3f 



0, 


Uj 






= 0, 


Uj 


= 1, 




0, 


Uj 


1 — ' 


^3 = 



(4.10) 



up to canonical transformations. The fourth solution ( IV) is the combination of the first and third 
solutions for the different j 's. 

Below we mark all the superintegrable systems by number of the solution from this Proposition. 
If we suppose that after point transformation 



X = 2:1(9), 

Px = wii((7)pi + wi2{q)p2, 



V = Z2{q). 

Py = W2l(g)pi + W22(q)p2, 



(4.11) 



the kinetic energy has some fixed form 

T = Y1 = Siii^: y)pl + gi2{x, y)pxPy + g22{x, y)pl, (4.12) 

then we will get some additional restrictions on the separated relations and the Stackel matrix. 
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For instance, if we suppose that 

then one wiU get the following algebraic equations 

W11W21 = (5^"^)^^ , W12W21 + W11W22 = 0, W12W22 = (4-13) 
and the partial differential equations 

{x,Px} = {y,Py} = I, {Px,y} ^ {Py,x} ^ {px,Py} ^0. (4.14) 

on parameters a and functions ^i, 2(91, 92), '^kjili, ^2)- 

The remaining free parameters 7,7j',7j' determine the corresponding potential part of the Hamil- 
tonian V{x,y). In fact, since integrals Hi,2 is defined up to the trivial shifts Hk Hk + Cfc, our 
potential V{x,y) depends on three arbitrary parameters only. 

4.1 Examples 

According to ^25j in this section we consider superintegrable systems on a complex Euclidean space 
i?2(C) with the following Hamiltonian 

Hi^PxPy + Vz{x,y), (4.15) 

where subscript Z = I, II, III or IV indicates on the type of the solution (|4.10p . The passage from 
the conformal coordinate system (x, y) to another coordinate systems can be realized by using the 
Beltrami partial differential equations. 

Remark 7 According to [211 HS] the Stackel transformations (|3.13p relate systems associated with the 
same solution of (|4.10p because the corresponding Stackel matrices Sz{q) are differed on the first row 
only. 

Example 1 Let us put ki = 1 and K2 = 1 in (14. 9p and try to solve equations (|4.13p - (|4.14p . Here is 
one system with first order additional integral Km (|4.5p 



vjn = li^y + 72(2; + y) + 73(3; - y), 
and seven systems with cubic integral of motion Km (|4.5p 



^(1) ^ _J]_ ^ 72 ^ 73 {y + x) 

T.(2) , 72 , 73 



[y-xY [y + xf 
(1) _ 71 , 72 , 732: 



^y{x-i) ^y{x + 1) VP 

T/(2) , + 1 73a: 
Vj\ = 71 a;?; + 721/ ' 



2 



Vx^+l V X 

y{2) ^ 71 72 ^ 73 

Vin = liy-'^^+l2xy-'/^ + j3X, 

T. , / ^ , laiSy - x){y - 3x) 

= (^T^ + ^^(^^") + 3 ' 

Solution of the recurrence relations (|4.7p 

Km-i = ±(2pip2UiU2 + 2viV2f + Ui52 + "251) (4-16) 

is the additional polynomial integral of motion, which is the first order polynomial for the system with 
potential vjjj and the second order polynomial in the momenta for another systems. 
All these systems are listed in the Drach papers [7j . 
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Example 2 At Ki = 1 and K2 = 2 there is one superintegrable system with cubic additional integral 
Km (|4.5p and quadratic integral Kg (|4.6p 

Viii = 71 (3a; + y){x + 3y) + 72(x + y) + 73(2; - 2/) , 

and seven systems with the real potentials 

72 , 73 



Vi = -fi{3x + y){x + 3y) + 

T.(l) , 72 ,73 

^11 -^^^y + ^^ + z^ 



{x + yY [x-yf 



_ 71 ,72 73y 

T^(3) _ 71 , 72 , 73 



'II 



xy x^/^y^/^ 



T/(4) , 722/ , 732/^ 

..(1) _ 71 , l2iVx-^) 73 



T/(2) I / \ , 73 

^iv = 712^2/ + l2{x-y)^ 



(a; + yY 

for which integrals of motion Km and Kg (|4.5M.6|) are fifth and sixth order polynomials in the momenta. 
Solution of the recurrence relations (14.71) 



Km-i = {Km, H2} - 4piP2MiU2 (2/1^2 + .92) + 4w2(2/t;i + .gi)(/w2 + 92) + (4//i2 + 52)^1 
= 4^2 (Ml ^2 +Pi)- (4//12 - 272' )Ai - 4/^231 . 

is the additional polynomial integral of motion, which is the second order polynomial for the system 
with potential Vm and the fourth order polynomial in the momenta for another systems. 

Solution of the next recurrence relation Km-2 = {Km-i, H2} is the rational function in momenta. 

Example 3 Now we present some superintegrable Stackel systems at ki = 1 and K2 = 3. Here is one 
system with cubic additional integral Km (|4.5p 



Vhi =Mx + 2y){2x + y)+-f2ix + 2y) + 73(22; + 2/) , ^ = ( ? J ) 
and seven systems with the real potentials 

72 , 73 



Vi = 7i(x + 22/)(2a; + 2/) + 



T/(i) , 72 , 73 

^11 = 712^2/ 



(a; + yY {x — yY 



Vii 



2-2/3^^4/3 ' j.l/3y5/3' 

(2) _ 71 _|_ 72 ^ 732/^ 

'xy x^/"^ x^ ' 

(3) _ '71 72 73 
xy a;4/3y2/3 + r^7/6y5/6 ■ 



Vii 

T/(4) , 722/^ , 732/^ 

Vjj =jixy+—^ + -^, 

Vifi = 71 ix' - 5x^ + 4j/) + ^ + ^ , 
Viv ^li{x + 5y){5x + y)+ j2{x - y) + 



{x - yY ' 



for which integrals of motion Km and Kg (|4.5M.6|) are seventh and eights order polynomials in the 
momenta. 
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4.2 Algebra of integrals of motion 

For all the superintegrable systems considered in the previous section the algebra of integrals of motion 
Hi 2 and Km, Kg (|4.5[) is the polynomial algebra in terms of the coefficients of the hyperelliptic curves 



{H2,Km} = 2Ki , {H2, Ki} = 2fKm , {Km, Kf} = ±$z , 
where polynomial Fz depends on the type of solution (j4.10p only: 

$7 - 2{^fh2 - glT'-^' (4/(k?/i2<?i + «2/»i52) - 5i52(«2.gi + 492)) , 

^11 = ^{^fh2-9lT^-^'[^f{n2 + ni)h2hi~Kihigl-K2h2gi)TK^m, 

$/// = A{A!h2-glT-'^'(^^f{nih2 + n2hi)-Kigl-K2gi)f , (4.17) 

$/y = 2{Afh2-glT--^'(^^f{2K2fhi-Kih2gi)-2^i2.f9l + ^igigi)- 

Here v ~ 1 for ki = K2 and v — 2 for ki 7^ K2. Choice of the sign + or — depends on k's too. 

The Stackel transformations p.l3p relate systems associated with one type of the solutions (|4.10p , 
whereas algebra of integrals of motion is invariant with respect to such transformations. 

5 Elliptic angle variables 

The first demonstration of the existence of an addition theorem for elliptic functions is due to Euler 
[8], who showed that the differential relation 

Ki dx K2dy 

■ =0 (5.1) 



connecting the most general quartic function of a variable x 

X ^a^^ + 46x3 g^^2 ^_ 4rfx + e (5.2) 

and the same function Y of another variable y, leads to an algebraical relation between x and y, X 
and Y at integer ki_2- 

As an example at ki_2 = 1 one gets 



x-y / 



a(x + y)-^ + 46(x + y) + C 



where C is the arbitrary constant of integration. This algebraic relation when rationalized leads to a 
symmetrical biquadratic form of x and y 

F(x, y) = axV^ + 26xy(x + y) + c{x^ + 4xy + y^) + 2o?(x + y) + e = 0, (5.3) 

which defines the conic section on the plane (x, y), which then will be classical trajectory of motion in 
the configurational space. 

According to [3] we could replace the constant C in the Euler integral relation by 4c + 4s, where 

f(x,y)-VTVF 1 ( Vx-yy V (x + y)^ 

s = TT, = T o(x + y) - c, (5.4) 

2(x-y)2 4\^ x-y y 4 

is the famous algebraic Euler integral. Treating s as a function of the independent variables x and y, 
one gets the following addition theorem 

dx dy ds ^ , ^, 
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Of course, the Euler addition theorem is a very special case of the Abel theorem [3]. As an example, 
the similar explicit addition formulas for genus two hyperelliptic curves may be found in [5] . 
The polynomial S in (|5.5p may be defined in the algebraic form 

_ (Fix + Y2)VX - (Xiy + X2] VY 

7 ^ , (5.6) 

where 

Xi = (ax^ + Sbx^ + 3cx + rf), X2 = 6x^ + Scx^ + 3rfx + e 
and similar to 1^1,2, or in the Weierstrass form 

5 = 4s3 _ g^s - 53 (5.7) 

where 

g2 = o,e — Ahd + 3c^, 53 — ace + 2bcd — ad^ — eJj^ — c^, (5.8) 

are the quadrivariant and cubicvariant of the quartic X [2] , respectively. 

For Ki = 1 and K2 = 2 we present Euler's equation of the trajectory of motion at c = d = only: 



F(x, y) = (e - ay4)x + 2y^e{ay^ + 6cy2 + e) , 

see page 355 in [5]. 

5.1 Classification of the Euler superintegrable systems 

According to [11], in order to use the Euler addition theorem (|5.5[) for construction of the superinte- 
grable Stackel systems we have to start with the genus one hyperelliptic curve 

^l^^P{\), where P(x) = X, 

and a pair of arbitrary substitutions 

^3 = ^3 il3 ) Aij = {qj )pj, j = 1 , 2 , 

where p and q are canonical variables {pj, qi} = 6ij. 

This hyperelliptic curve and substitutions give us a pair of the separated relations 

p2 y2(^^.) ^ avjiqjf + ibvjiqjf + GcVjiqjf + ^dv,{q,) + e, 1, 2, (5.9) 

where coefficients a, 6, c, d and e of the quartic X (|5.2p are linear functions of integrals of motion Hi 2- 

a = aiHi + a2i?2 + a, 6 = + P2H2 + /?, . . . , e = eii?i + £2-^2 + e. 



As above the separated relations (|5.9p coincide with the Jacobi relations for the uniform Stackel systems 

[Miiniiii] 



P3 



Y.H,Sk,-U,iq,) (5.10) 

\ k=l 



where S is the so-called Stackel matrix and Uj is the Stackel potential: 

Sij = uJ^{a^vj + A(3ivl + + ASiVj + e^), 

«,j = l,2. 

Uj = Uj ^{avj + 4l3vf + 6jv] + 4Svj + e). 

Solving these separated relations (|5.9ll3.3p with respect to Hi_2 one gets pair of the Stackel integrals 
of motion in the involution 

2 

Hk = J2iS-^)jk{p^j+U,iq,)) , A:- 1,2, (5.11) 

3 = i 
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and angle variables 

canonically conjugated with the action variables Hi 2- 

In generic case the action variables (15.12^ are the sum of the multi-valued functions. However, if 
we are able to apply some addition theorem to the calculation of uj2 

1 S2i(A)dA 1 r^(«=) S22(A)dA 1 r As 

UJ2 = ■ ' 



/P(A) 2j ^/p{Xj 2j 75' 

then one could get additional single-valued integrals of motion s and S: 

{Hi,LU2}^0 ^ {Hi,s} = {Hi,S} = 0. 

Of course, opportunity to apply some addition theorem to computation of the angle variable lu2 and 
of the additional single-valued integrals of motion s, S leads to some restrictions on our quartic P(A) 
and substitutions A^ = Vj{qj), fij = Uj{qj)pj, see [TT] . 

For instance, if we want to use the Euler addition theorem (|5.5p we have to put 

S2l(A) = l, S22(A)=±1. 

Such as 



Uj{qj) 

these restrictions are equivalent to the following equations 



' - a2w| + 4/32w| + 672w| + 4(52Wj + £2 , ki = 1, ^2 = ±1, (5.13) 



J — LX2 '^j ~ 



on functions u{q),v{q) and coefficients a2, P2, ■ ■ ■ , £2 of the quartic, because we have to solve these 
equations in some fixed functional space [S^ . 

It's easy to prove that there are five monomial solutions 



I 


Uj 


= 1, 




= 9i> 


£2 y^O 


II 


Uj 


= qj^ 






(52^0, 


III 


Uj 


= 1, 




= q-\ 


72^0 


IV 


Uj 








/32^0, 


V 


Uj 


-qj^ 




-qj\ 


a2 7^ 



(5.14) 



up to canonical transformations of the separated variables {pj,qj) and transformations of integrals of 
motion Hj — ^ crHj + p. Here notations q;2 7^ 0, /?2 7^ 0, . . . mean that other parameters are equal to 
zero. 

As above, if we suppose that after some point transformation (|4.1ip kinetic part of the Hamilton 
function has a special form ()4.12p . then one gets some additional restrictions ()4.13M.14| ) on the 
coefficients of the quartic P(A). 

The remaining free parameters a, . . . , e determine potential part of the Hamiltonian V{x,y). In 
fact, since integrals Hi 2 is defined up to the trivial shifts Hk — > + Pk, our potential V{x, y) depends 
on three arbitrary parameters only. 

Summing up, we have proved that classification of the Euler superintegrable systems on the plain 
is equivalent to solution of the equations (|5. 1314. 1314. 14p . For all the possible solutions additional 
quadratic in momenta integral of motion looks like 

F{vi,V2) - P1U1P2U2 , 

K2 = s + c= — h c, 

2[vi - V2Y 
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see ()5.4p and additional cubic in momenta integral of motion is equal to 
see (EJI and (l57l). 



Remark 8 Of course, additional integrals of motion for the n-dimensional superintegrable systems are 
related with the Abel theorem as well. As an example we present the quadratic in momenta Richelot 
integral [Hj 

for the uniform Stackel systems linearized on the Jacobian of the following hyperelliptic curve 

a2riA^" + h fliA + flo, 

such that the Stackel matrix is a lowest block of the corresponding Brill-Noether matrix [201 HI] ■ Here 
F(A) — Y\{\ — \j) and \j = Vj{qj), /ij ~ Uj{qj)pj are the suitable substitutions. 

5.2 Examples 

Let us find all the Euler superintegrable systems on a complex Euclidean space E2 (C) 

Hi=p^Py + V{x,y) (5.15) 

with the real potentials V . Solving equations (|5. 1314. 1314. 14[) one gets the following five superintegrable 
potentials 

Vi = a(a; + y)+/3(y + 3a;)x"i/2+72;~^/^ 
1^2 = ay(a; + y2) + /3(a; + 3y2)_^^y^ 

^3 = a(a; + 3y)(3a; + y)+/3(x + y)+ ^ 



(x - yf- 
V4 = axy~^ - (3y~^ - -fxy, 

^ x-^ ~^ x^^^y/y — 1 ~^ x^/'^y/y — 1 ^ ^ 

Recall, that implicitly all these systems have been found by Euler in 1761 [8]. Potentials Vi and V3 
in explicit form have been found by Drach 7], the {£) and (g) cases, whereas potential V2, V4 and V5 
may be found in [14^ , as the Eiq, and £'17 cases, respectively. 

The first and fifth solutions (|5.14p of the equations (|5.13p are related with potentials Vi 2- The 
second and fourth solutions give us potential V3, whereas third solution yields potentials V4.5. For 
brevity we present the main stages of calculations only, all the details may be found in [llj . 

Case 1 If we take first solution from (j5.14p 

= ±1, Vj = ±qj, => p\ = P{qi), {~p2f = P{~q2) 

and quartic 

P{\) = -|a4 + + Hi\^ + 27A + H2, 
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a{x-yf P{x~yf l{x - y) 



then after the following change of variables 

('71-92)^ {qi + q2f P1+P2 
x = -. , Px ^ Pi - P2qi - q2, y = 7^ > Py — : — ^- (5-17) 

4 Z qi + qz 

we obtain the Stackel integrals 

/ N Pi^x + y) 1 

Hi = PxPy + a[x + y) H -j= h —j= , 

•y/X \/X 
H2 = [Px - PyjiPxX - Pyy) - 

the Euler integrals of motion 

2 

P,, ax ^ 

J. /t; P\kPx-Py) a{p^-py)x f]{2p^x - 3pyX + Pyy) jpy 
4 + 2 + + 

The same system may be obtained by using fifth substitution from (j5.14l) . 

Case 2 Using the same first solution (|5.14p and another quartic 

P(A) = -jX^ - I3X^ - + HiX + H2 

after the following change of coordinates 

(91+92)^ P1+P2 91-92 

x = , Px^ ■ , y= , Py=Pi-P2, 

4 qi+q2 2 

we can get superintegrable Stackel system 

Hi = PxPy + ay{x + y'^) + (3{x + 3y'^) +'^y, 

„2 



Py ^ 2 ^ a{3y^ + x){x-y^ 2.^1, 
-f+PxX- yPyPx + ^ + 2f3y{x -y ) + -( 



with the quadratic Euler integral 

4 4 2 

and with the cubic integral of motion 

pI a{3pAx + y')~Pyy) p^Qp^y - py) 
K, = VS=- + + + — . 

The same system may be obtained by using fifth substitution from (|5.14p as well. 
Case 3 If we take the second solution from (|5.14[) 

"j=9j> '"j^'ij' ^ plql^P{ql), plql^P{ql) 

and quartic 

P(A) = -aA" - |a3 + HiX^ + H2X + 7, 
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after canonical transformations (I5.17P we will get the following superintegrable Stackel system 

Hi = PxPy + a{x + 3y)(3a; + y) + (i{x + y) + 



{x - yf ' 

H2 = (Px - PyKPxX - Pyy) - 2a{x + y)[x ~ y) 

2 (x - yY 

Using the Euler addition theorem ()5.5|) we can get the Euler integral 

K,^s + c=^P^±^+a{x + yr + ^ix + y) 
lb 4 

and the cubic in momenta integral of motion 

iPx-Py)^{px+Py) ~ 2/)(p^(5a; + 'Sy) - [ix + 5y)py^ 



32 

P{Px~'Py){x~y) liPx+Py) 



16 ^{x-yY 
The same system may be obtained by using fourth substitution from (j5.14p . 
Case 4 If we take the third solution from (|5.14p 

and quartic 

P(A) = aX^ - f3X^ + H2\^ + HiX - 7, 
then after canonical transformation 

Pigi + 92^2 _ gi - 92 _ V^i92(pigi - 92^2) 



X = y^qiq2, Px = ; — , V = , Py 



V^iQ2 Vm2 gi + q2 

we will get the following superintegrable Stackel system 



rj , "2/ , /? , rr 2 , (PxX - yPyf a{y^ + 1) (3y 2 
Hi = PxPy + — + — + ixy, H2 = Py H — ^ +JX . 

Using the Euler addition theorem (|5.5p one gets additional quadratic in momenta Euler integr 

„ ^ (PxX - ypyY ay'^ Py 

^ * 16 4a;2 4a; 

and the cubic integral of motion 

- pliPxX-yPy) a{3ypy + xp^) Ppy j{pa;X - ypy)x'^ 



8 4a; 8 

Case 5 Using the same third solution (|5.14p and another quartic 

P(A) = 4aA* + 4/3A^ + H2}? + 27A + i?i 

after the following change of variables 

_ 9192 _ Pigi + 92^2 _ 9i_+9| _ gig2(pigi - g2P2) 

X — , Px I?/ „ - Py 2 2 ' 

we can get more complicated superintegrable system with the Stackel integrals of motion 

rr _ a (3 7 

rr ( V ^ A ^ay 2/?(2y-l) 27^^^ 

H2 = {XP, - Py - Pyy){xp^ + Py - Pyy) - — ~ -^JJ^== - ^==. 
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The Euler integral of motion is equal to 



A2 = S + C = -i= , 

whereas the cubic in momenta integral reads as 

Py{xPx + Py - VyVf , '^Oi{y ~ l)py f3{p^x - 3pyy + 3py) -fx^/'^{pxX + py - pyy) 



This list of examples may be easily broadened because Euler proposed construction of the algebraic 
integrals of motion for the equations (|5.ip with any integer k's. We discuss the Euler construction 
of the algebraic integrals of motion for superintegrablc systems at ki,2 = ±1 only. Of course, the 
most interesting case is the case with the different ki_2 in (I5.1|l for which the corresponding additional 
integrals of motion will be higher order polynomials in momenta as for logarithmic angle variables. 

5.3 The quadratic integrals of motion 

It is easy to prove that the algebra of integrals of motion Hi,2 and K2 is the quadratic Poisson algebra 
because 

Vi, 



{H2,K2} = aK3 = a^S{s), where { a = -2, 1/2,^4,^5, 

^3, 




and 



2 2 

{H2{H2,K2}} ^ {H2,<J^} ^ ^ S' ^ ^[V2s'' ~ g^) = i/2, ^^2), (5.19) 

where H2, K2) is the second order polynomial such as s = K2 — c and 52 is quadrivariant of the 

quartic (|5.8[) . Another details on the quadratic Poisson algebras of integrals of motion may be found 
inia. 

The search of the two dimensional manifolds whose the geodesies are curves which possess two 
additional quadratic integrals of motion was initiated by Darboux [i], who found five classes of the 
metrics. These metrics or "formes essentielles" are tabulated in "Tableau" by Koenigs [16] and in p3] . 

The superintegrablc Darboux-Koenigs systems have a generic conformal Hamiltonian 

where the Darboux-Koenigs metric g is a metric on the Liouville surface [H [S] if 

g(C,^)=i^(e + ry)+G(e-ry), and K2 = p^"" + p,,^ - 2p^p^^^ + Q{i,^) 
or metric g is a metric on the Lie surface if 



g(^,^)=^J^W+G(77), and K2 = p^^ - 2p^pr,^%^ + Q{i,Ti). 

Superintegrablc systems associated with the Liouville surfaces are separable in the two different or- 
thogonal systems of coordinates. It means that two pairs of integrals of motion {Hi, H2) and {Hi, K2) 
take on the Stackel form (|3.2p after some different point transformations (|4.11[) . 

For the superintegrablc systems associated with the Lie surfaces only one pair of integrals H2) 
may be reduced to the Stackel form (|3.2p . whereas second pair of integrals {Hi, K2) doesn't separable 
in the class of the point transformations (14. lip . 

It's easy to prove that two systems with potentials Vi and V2 are defined on the Lie surfaces. The 
remaining systems with potentials ¥^,¥4 and V5 are defined on the Liouville surfaces. The second 
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separated variables (71,2 for integrals of motion Hi,K2 may be found by using the software proposed 
in [in]: 



qf-qi _qi + qi 



for V3, 




x = q2qi^, y = q2q\, forV4,V5. 

It is easy to prove that the corresponding separated relations p| = P{qj) define two different zero-genus 
hyperelliptic curves and lead to logarithmic angle variables ^Sl- Thus, for these three Eulcr super- 
integrable systems there are two different addition theorems: addition theorem for elliptic functions 
(|5.5|) and addition theorem for logarithms Inx + Iny = Inxy. So, multiseparability of the supcrin- 
tegrable systems may be associated with occurrence of the different addition theorems for a given 
superintegrable hamiltonian. 

According to [13] , the Stackel integrals of motion are in the bi-involution with respect to a pair of 
compatible Poisson brackets 

where the second Poisson tensor P' = NP is obtained by the following recursion operator N: 

Opk ^ Opi 

Here L is so-called Benenti tensor [T]. There is symbolic software [T^, which allows us to obtain this 
tensor starting with the Hamiltonian Hi only. 

So, it is easy to prove that any superintegrable system associated with the Liouville surface has 
two different linear in momenta Poisson bivectors P' and P", which are compatible with the canonical 
bivector P and incompatible to each other. On the other hand, superintegrable system associated 
with the Lie surface has one linear Poisson bivector P' only. It will be interesting to understand the 
geometric origin of this phenomena. 

Moreover, we can find some integrable systems associated with the Lie surfaces, which doesn't 
separable in the class of the point transformations and which hasn't linear second Poisson bivectors. 
As an example, we present the following integrable system with the quadratic integrals of motion 

Hi = p^py + ay-^^^ + (3xy-^/^ + jx, K2 = pI - Ajiy^/^ + 2jy , 

which hasn't any other quadratic integrals. It will be interesting to classify such systems and describe 
the non-linear bi-hamiltonian structures associated with the Lie surfaces. 
The research was partially supported by the RFBR grant 06-01-00140. 
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